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FOREWORD 


This report covers the second portion of the notes prepared for a seminar, "Mechanics of 
Composite Materials, "presented at the Air Force Materials Laboratory in April and May 1966. 
v The work was initiated under Project No. 73401 "'Nonmetallic and Composite Materials," 
Task 734003, "Structural Plastics and Composites." The seminar consisted of Part I - 
Introduction, and Part II - Mathematical Theory. 

The manuscript of this report was released by the author in June 1966 for publication as an 
RTD Technical Report. 

■> • 

This technical report has been reviewed and is approved. 

7T 5 

R. T. SCHWARTZ, Chief*" 

Nonmetallic Materials Division 
Air Force Materials Laboratory 
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ABSTRACT 


This report covers some of the principles of the mechanics relevant to the description c 
composite materials. The contents of these notes may provide useful information for tk 
understanding of current publications and reports related to composite materials. 

i-- Emphasis is placed on the use of lndicial notation and operations. The rules governing th 
use of the contracted notation are also outlined. The generalized Hooke’s law and its train' 
formation properties, material symmetries, and engineering constants are also discusser 
The plane strain and plane stress problems are discussed in detail. Finally, the elastic modu> 
of laminated anisotropic materials, and the strength of both unidirectional and laminated com 
posites are covered. 
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Each tensor can also be arranged in a matrix form. The indices of each teneorial compo¬ 
nent can be associated with a specific position in the matrix. For example, a vector can be 
expressed by A, or a column or row matrix form as follows: 






where i * 1, 2, 3, which indicates that the space is three-dimensional. In 2-space,= 1, 2, 
Equation 2 becomes ' ' 



(3) 


For a second-rank tensor a 
are, respectively: 




the corresponding matrix forms for 2-space and 3-space 



For a fourth-rank tensor, the matrix form In 2-space is: 




0 . " 

o 

0 . 

0 „ 



1 III 

list 

Ill* 

118 1 


Vi* 

a ttn 

*8888 

*8818 

°888 I 

6 . ■ •' 

*•811 

*1 888 

*1818 

*1881 

'■ . 0 


.*8111 

*8188 

*8118 

*8181 


A fourth-rank tensor in 3-space contains 81 components which can be arranged in a 9 x 9 
matrix. c- 


Matricea contain an array of numbers. The numbers arranged in appropriate positions in 
a matrix may represent the components of a tensor. This does not mean that tensors and 
matrices are identically equal. The components of a matrix may be arbitrary and completely 
unrelated. But the components of a tensor, whether in the indidal notation or the matrix form, 
are not arbitrary. The components are governed by a set of rules, called the transformation 
equations. Therefore, a tensor can be considered as a special type of matrix. Tensors are 
in fact defined by the transformation equations. If a set of lumbers satisfy certain transfor¬ 
mation equations they are by definition components of a tensor. For each tensorial rank, there 
is a corresponding transformation equation. Thus, the coordinate transformation is a basic 
feature of tensors. 
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In a simple geometric term, a coordinate transformation Involves the rotation of one 
reference coordinate system relative to another. The relation between the old and new co¬ 
ordinates, designated by x { and xj, are, for 1 * 1, 2, and 3 as follows: 


*1 S *.» *. + ».* *2 + *13 *3 
*2 * *21 *1 + *22 *2 + *23*3 


*3 * V *1 + *32*2 T *33*3 


where t^ are direction cosines between the 1-axis in the hew system and the J-axis of the old 
one. Equation 7 may be written in a matrix form: 


r .. 

*12 

*13 

f s. 

*22 

*23 

31 

*32 

*33 


where the usual rule of matrix multiplication applies and this rule can be expressed in the 
indtclal notation as follows: 


*1 * *11 *1 * »». *1 + *.2*2 + *.3 ' 


1*1 

* • I 

* i*i 


*2 * ui *2* “i * *1 + ^2*2 + ^23*3 


*3 * Si *1 * 'si *1 + * 31*1 + f 33*3 


T’vo useful conventions of the indictal notation can now be introduced: 

V . . ■ . 

1) Range Convention: v " '■> 

Unrepeated index (free lnd» '.) takes all the values 1, 2,.... n, where n is the dimension of 
the space. \> 

2) Summation Convention: w ' . .. 

Repeated index (dummy index) calls for the summation with respect to that index with the 
range of summation determined by the range convention. 

With these conventions, Equation 9 can be written as: 


a. * t. 

I M J 
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j where i is the- free index and is equal to 1, 2, and 3 in accordance with the range convention, 

J J is. the dummy index, as it is repeated in the right-hand side of Equation 10, and a sum- 

c mation of j with the same range as i, i.e, 1, 2, and 3, is implied. The use of elimination 

j convention replaces the summation sign (capital Sigma) in Equation 9. Both range and sum- 

£ mation conventions are applicable to all square matrices (because all indices must have the 

I same range) and are therefore not limited to tensors. 

.vEquations 7, 8, and 10 are identical equations and it is quite clear that the use of the in- 
dicial notation introduces a significant simplif. cation. 

TRANSFORMATION MATRIX 

As stated earlier, t^, which for convenience is designated T, are direction cosines between 

the old and new coordinate systems. Knowing the angle of rotation between the two coordinate 
systems, the direction cosines can be determined immediately. The components of the trans¬ 
formation matrix T for a rotation about the 3-axis can be derived from Figure 1 and are listed 
in Table 1: ~ 


w 



I I 


Figure 1. Coordinate Rotation 
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TABLE 1 

1 -MATRIX FOR ROTATION ABOUT 3-AXIS 


Components 

Angles 

Direction Cosines 

*11 

ct> 

u 

fH 

o 

"h 

cos 9 

*12 

1* 02 = <f) 

sin 9 

*1.3 

1* 03 = 90° 

0 

*21 

2* 01 » 90* + 9 

-sin 9 

*22 

: 2 02 = 9 

cos 9 

*23 

i 03 = 90 s 

0 

*31 

3'01 * 90* 

0 

*32 

3‘02 = 90* 

0 

*33 

3*03 = 0 

1 


If m * cos 9 and n * sin 9, the T-matrix above can be arranged In a matrix form as follows: 



T « t 



m n 0 ' 
-a m 0 
0 0 I . 


(ID 


For a rotation in the opposite direction at that shown in Flgare 1, 9 is negative, and the 
transformation matrix, designated T", is as follows: 


T“ » 


m -n pi 
s m 0 


( 12 ) 


0 0 I 


The components of the T-matrix for any other rotation can be derived similarly from 
angles listed in Table I. The subscripts of the t^ refer to the new and old axes, i.e., 

angle for the direction cosine is that between the i-th axis of the new coordinates and 
j-th axis of the old coordinates. 

There are a number of features of the T-matrix, as follows: 


1) When 9-0. 



( 13 ) 
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f , , 

Tits is called the identity transformation, for which x! = x.. The new and the old coordinates 
are identical. , . , 

2) When 6 * 180 V : , 

I - "' 0 °1 

T = I 0 -1 0 I (14) 

L 0 0 I J 


This is called the central inversion, for which xj = -x„, x^ - -x 2 , and xj = x<j. 

3) A proper rotation is defined by a T-matrix with its determinant equal to +1, i.e., |t| =+l. 
For a 3 x 3 determinant, its numerical value is (he sum of the following products: 


n 



f .l 

*22 

*33 

♦ 

♦ 

t 

12 

23 

31 

t 

t 

♦ 

21 

32 

13 

t 

t 

t 

31 

22 

13 

t 

t 

t 

21 

12 

33 

t 

♦ 

t 

1 1 

23 

32 


(15) 


The determinants of the T-matrix in Equations 11, 12, 13, and 14 are equal to +1, and are 
thus proper rotations. Geometrically speaking, all of these transformations preserved the 
right-hand system of coordinates, as shown in Figure 1. This system satisfies the right-hand 
rule, when it is applied to the coordinate axes, by rotating the 1-axis toward the 2-axis with 
the 3-axis as the advancing screw. 

■ ■ i 

If a transformation changes the right-hand system of coordinates into a left-hand system, 
this is called an improper rotation, for which |T| ■ -1. An example of this transformation is 
the reflection of the 1-axls, such that xj - -x^, xj * x , and xj » x^. The T-matrix becomes: 



Note that | T| ■ -1, which Indicates that the rotation is improper and the new coordinates x! 
are now a left-hand system. This is shown in Figure 2. rr 







The Kronecker delta In a matrix form is: 

[ I 0 0 1 

0 10 (19) 

0 0 ) J 

which is also a unit matrix. The use of Kronecker delta is a very important tool in tensor 
operations. 

TENSOR CALCULUS j ‘ ! 

4 

Three salient features of the indicia! notation mentioned thus far are the range convention, 
the summation convention and the Kronecker delta. When they are applied simultaneously, one 
can show that a 



.t 
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8.. = 2 fori * I, 2, 

ii 


* 3 for i * I, 2, 3 


( 20 ) 


Another feature of the Indicl&l notation is a direct correspondence between the number of 
free indices, k (not the dummy index, because a^ is a scalar, a^, a vector), and the tensorial 

rank. The number of components, N, in each term is determined by Equation 1, where k is 
the number of free indices and n is the number of dimensions. 


Finally, in this notation, a comma represents differentiation, as follows: 

da. 

IT s Vr 

1 


( 21 ) 


By using tiie summation convention, one obtains 

do. 


a i *1 * dx 


dOj do 2 do 3 

dx dx dx 
I 2 3 

» V* a 

* dlv o 


( 22 ) 


where, in the last two steps, the oooventiooal vector notation at divergence is used. Similarly, 
for a scalar function A, 


A,; 


« 

dx 

I 

, JA. +. -JiL. 

da, d« 2 

* VA 
■ grad A 


iA_ 

dx- 


(23) 


where, in the last two steps, the conventional notation for gradient Is used. 
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Finally, 



d 2 A 

d«A 

d 2 A 


d* 


d Z A 
dx 2 


dlv ( VA) 

V 2 a 


a 2 a 
a*! 



which is the Laplacian operator. Since the tensorial rank can be determined by the number of 
free indices, and A,^ are scalars; and A^, a vector. The fact that the tensorial rank can 

be determined by observation is a feature of the indlcial notation that does not exist in the 
conventional vector notation. 

The divergence theorem can be written as follows: 

/ a.n.dA * / a. ,dV (25) 

* A ' ' V 


where V = volume, A » surface,n t • exterior normal to A. Substituting a^ * A tJ into Equation 25 

/ A, .n.d A * / A . ,dV (26) 

A ' V ' 


This relation can be used, for example, in the derivation of the Fourier heat conduction equation. 
Equation 25 can also be generalised to a vector equation by using in place of a^ 

/«r t| n.dA • /<r. ,dV (27) 

A 1* 1 y l 

This relation will be used In the derivation of the equilibrium equation. 

STRESS TENSOR * .. . 

Stress is a measure of the internal foroes in a continuous medium induced by surface forces 
applied to. a body. The relationship can be defined by 

T. « tr. ,n. (28) 

i j» 1 


where T { = surface traction, * stress tensor, and n { = directional cosines of the surface 
on which Tj acts. Equation 28 in 2-space is: 


T , '‘'ll". + ‘Vz 


T 2 "kV* ff 22 l *2 


(29) 
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The relations can be seen in Figure 3. 

2 



Figure 3. Stress Tensor in 2-Space 


The normal component of the surface traction, T q , can be obtained by 

T » ■ V. ’ »|.V. 


(30) 


In 2-space ^ 

T n ««r„» 2 | f 2,r l2 B ."2 +<r 22 B2 2 

The tangential component of the surface traction, T ( , can be obtained as follows: From the 
directional cosines n } of the normal to the surface, one can find the directional cosines t ( in 
the tangential direction by transforming n^ through 90*, Since n t is a first rank tensor. 


(31) 


where 


* Vi 



n 

m 

0 


0] [0-1 o' 
0 * 10 0 
I J L 0 0 I . 


(32) 


(33) 
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Thus ■ 

V'W n I ' *3 * n 3 (34) 


Hence 


T, * T. t . = or. n t. 
t i i ji j • 


(35) 


Fur 2-space 


T. t a n t + n I- + <r n t + oint 
* II I I 12 I 2 21 2 I 22 2 2 

s-ff n n 4ff n 2 — & n^+crnn 

1112 12 I 21 2 22 2 I 


(36) 


Since 


T , ’ T n ”, - T , "2 

T 2*~ T n * l + V* 


(37) 


From Equations 31 and 36 


T l * V* + **i2*iS , 


(38) 


Similarly, 


T 2 1 V* "2 + *'l2"l"2 2 + '22*2* ~ *2 " 

+ °22 B | 2 "2 * "’.2-. *'V* 


(39) 


Equations 38 and 39 agree with Equation 29, as expected. Finally, one can easily show that 


+ T , 2 * <v/ ♦ < T .v* 


"V, M i V !,IT ,' l tT ! l 2' ! 

H.n.t 


I I 

’ T , 2 » T 2 2 

. T* 


(40) 
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The equilibrium equation of a continuous medium can be derived from the conservation 
principle of linear momentum, as follows: 


/ T.dA + / B.dV « fp u'.dV (41) 

A * V 1 V 


where = surface traction, B } = body force, p = density, u t = displacement, double dots = 
occoud derivative with respect to time, A = surface area, V = volume enclosed by A. 

From Equations 27 and 28 

/ T.dA * /cr.n.dA * / <r dV (42) 

A 1 A »• » V )'•) 

Since the volume is arbitrary, 

<r.. . + B, * nil. (43) 

j'.j i r • 

In the absence of body forces and time effect, the static equilibrium equation is 

or * 0 (44) 

J'i I 

From the conservation principle of angular momentum that if the resultant moment due to 
body and surface forces vanish and there is no other internal or external sources of moment, 
it can be shown that the stress tensor must be symmetric, i.e., 

or. . * <r.. (45) 

•I 11 


In view of the symmetry property of stress, Equation 44 can be written as 



This equation in 3-space can be expanded, as follows: 


**ll 



d<r \i 

d *l 



d 's 

£fl2 

.A, 


d '\ 

▼ 

d *2 ^*3 

is* 



a<r 33. 

da, 



^*3 


(46) 


(47) 


STRAIN TENSOR 

Strain is a purely geometric relation. It is Independent of the constitution of a body as long 
as it remains continuous during the process of deformation. 


% 
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Let points P and Q, separated ty a distance ds in an undeformed body, displaced to P 1 and 
Q 1 during a deformation process. The relations between,these four points are shown in Fig¬ 
ure 4. ■: 



0*( «' s + dx’j) 


Figure 4. Displacements of Points During Deformation 


From Figure 4, 


u i “l " *1 


(da) 2 * 8,jdX|dij 
(da*) 2 * S^dx’jda'j 


(48) 

(49) 

(50) 


If the process of deformation is referred to the final or deformed state (the Eulerian co¬ 
ordinates) 


Then 


Strain tensor, e^, is defined by 


.2 


« jiU*,) 


(51) 

1 Yi*"'} 


(52) 

WW^'i d *j 


(53) 


• 




(54) 

) 2 * 2«.,dx'.d*‘, 


(55) 
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By substituting Equations 53 and 54 into 55, one obtains: 


Since 


2s,. • 8.. - x_ ,x_ . 
i j i j m,i m.j 


vax - u 
*m, m m 

x - 8 . — u 

m t i u mi ''m,! 

X_ j - 8- I - u_ J 

m,j mj m,j 


(56) 


Then 


x •'*— 5 ( 8 . - u .) (8 . - u ,) 
fn,i mj mi m t i mj mj 

* 8. . - u. , - u. ,+• u_ . u_ . 

ii I,i j,l m,l m,i 


(57) 


By using Equation 56 


2s,. « u. . ♦ u. , -u_.y. . 
ij i,j J, • si,) 


(58) 


The Eulerian strain when expanded becomes 


du 

*11 * d%‘ ~ T 

I 



/ du_ 

d u_ \ 2 

{ _ 2\ . | 

_ I\ 

v I * 




(59) 


and 


± / 4 “. 4u 2 \ JL/""' -Hi* 

12 2 \7V377r 2 \TT t »v, TI 


d u. d u. du^ du 2 

1 - r-r 

I 4 *2 


d u_ d u 


U 3 gU 5 \ 

*'z> 


(60) 


In a similar way, one can derive the Lagrangian strain which refers to the initial or un¬ 
deformed state, as follows: > 


2s,. » u, , + u, . + u. u_ . 

I j I#) Jt 1 J 


(61) 


In expanded form 


d U| 

+ — 


SU 2 f . / 

4 “> f ■ 

dx, 

+ 2 

_l/ + \ 

d X , / ' 

a«, / _ 


(62) 
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and 


12 




d u. d u. 


3 u 2 * l 2 


aU 3 du 3 \ 
\ d *2 ' 


i» 2 


dx ( dn 


(63) 


The physical significance of the strain components can be illustrated by two examples: 
1) Normal strain. 

Let 

da* * dx‘ ( , implying dx’^ 5 dx* 3 1 0 


ds * (I — E ( Ida* * tl~E | >d*'' | 


(64) 


where 


_ ds - d» 

« 5 ds’ 


(65! 


From Equation 55, 


(d » 2 - (I- E l*'(.da* I* * 2# ,, d *‘, 2 


( 66 ) 


therefore 


•h ‘ E . -*i '■* 


(67) 


Thus, for Infinitesimal displacement, 




1661 


which means that, in the case of Infinitesimal strain, the normal component of strain cor¬ 
responds to axial elongation. For finite atrain. Equation 67 will apply. 

2) Shear Strain 

Let ,, . 

da* « da’j , with dx’ ( «di‘j*0 (69) 

da’ « da'j , with da'j * dx’^ « 0 (70) 


15 
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The angle 8 between the two vectors is 


From Equation 56 


Let 


cos 8 - dx.dx, / dsdi 

■ i 

= x l,2*i,3 d *‘2 d *3 /ds d7 

- jt 

cos 8 - (^S^ — ^*23*^*2 *** 3 


4> - 9 - 90* 



(71) 


and it represents the change in the angle between the two elements which was 8 in the initial 
state and became perpendicular in the final state. 

" f rom Equation 71 


sin da' 2 d*' 3 / ds ds 


(72) 


From Equation 65 


sin ^ * 2« 23 ✓ (I - E 2 > {i - E 3 » 


(73) 


F rom Equation 66 

•in ^ * 2« 23 /yil - 2« 22 )(l - (74) 

For infinitesimal strain, 

> $ « 2« 23 (75) 

Thus, the infinitesimal shear strain corresponds to one half of the change in angle from 90 s . 

r, ’ \ 

The strain-displacement equations for both Eulerlan and Lagrangian strains reduce to the 
same equation in case of infinitesimal strain as follows: 

#. ,* -f <«. . + «. J (76) 

* I 2 • ,| j,i 

This strain is a second-rank tensor. The conventional engineering strain is different from the 
tensoria) strain in the shear strain components. The engineering shear strain is related to the 
total change in an angle, rather than one half of the angle according to Equation 75, For this 
reason, engineering shear strain is twice the tensorial shear strain. 


16 
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SECTION II 

TRANSFORMATION EQUATIONS 


Tensors are defined by the following transformation equation: 


C ,, ... s I. t. t *, C 

im jn... ko Ip mn...op 


(77) 


where C....... is the transformed C ... : t.., the transformation matrix. In first, second, 

ij kl mn op i] 

and fourth-rank tensors. Equation 77 can be specialized as follows: 


C. = t.C'- 

i im m 

C'. . = f. t. C 
i j in jn mn 

C* - t t I t C 
ijkl im jn ko Ip mnop 


(78) 

(79) 

(80) 


FIRST-RANK TENSORS 

Equation 78 represents the following simultaneous equations in 2-space and 3-space: 


(81) 



r c ', n 


" *. i 

*1.2 


' C 





Z 







c 1 


t 

t 


C 



L 

2- 


21 

22 J 

l .. 


. 2 J 

r 

a 


" 



- 



c 

1 


*1 1 

*12 *13 


C . 


X 

*21 

*22 *23 


C 2 

c 

* 


i 

» f 



C 

L 

3. 


L 

31 

32 33 . 


L 3 


(82) 


Assuming that the transformation is a positive rotation T , which has the components shown 
in Equation 11, Equations 81 and 82 become 


C 

L 2 J L 


m 

- n 


n 

m 




< 


2 


. C 3j 



m 

- n 


’I 


(83) 


n 0 

m 0 

0 I 


C 2 

J 

1 

’«! " 


C 2 


1- 

o 

w 

L_ 


(84) 
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These are the transformation equations of T + for the first-rank tensors (vectors). 
SECOND-RANK TENSORS 

For second-rank tensors, the transformation Equation 79 in 2-space can be expanded as 
follows: 


1) When i = 1, j = 1, 


Summing m gives 


^ 11 " *lm , l n ^ mn 


C st t C. + t „ t C 
II l I I n I n I 2 In 2 n 


(85) 


Summing n gives 


C ,. S + # ,2 C ,2 » ^ , .2 (, .l C 2. - *.2 C 22» 


2) When i = 2. J = 2, 


C 22 * , 2m , 2n C mn 


Summing m gives 


C* « C. + t, C_ 

22 21 2n In 22 2n 2n 


. 86 ) 


Summing n gives 


C 22 * *21 tt 2» C l I + *22 C I2 * * , 22 (, 2l C 2I + *22 C 22 ’ 


3). When i » 1, J ■ 2, 


Summing m gives 


®I2 * , l«i*2« C mn 


C' • i M I, C, ♦»..»• C % 
12 II 2s In 12 2n 2n 


(87) 


Summing n gives 


C* * t (t C-’ *tC) + t(t C 4-10 ) 

12 II 1 21 l I 22 12 12 21 21 22 22 
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4) When i = 2, j = 1, 


Summing m gives 


Summing n gives 


C 21 ” *2m*ln ^mn 


C„, * i, C, 4 t„ t, C„ 
21 21 In In 22 In 2n 


C Z, ’ <1, "u C u + ',2°,2 '* , 22 ,, M C *I + , , 2 ‘W 


From Equation 11, the components of T-matrix in 2-space are 

T 


•(" •) 


By substituting into Equations 85 through 88, one obtains 


cl, 

« m 2 C, . + n 2 C_, 

♦ mnC, * 

4mnC-, 

11 

II 22 

12 

21 

c' 

* n 2 C +ml 2 C-- 

-mnC . 

- mnC , 

22 

N 

M 

12 

21 

C 42 

« -mnCj ( 4 mnC 22 

♦A* 

*•%« 

C 2I 

■ -mnC |( 4mnC 22 

- A,j 

♦A*. 


In a matrix form Equation 90 becomes 


V 


c « 

' 0 

C 'l2 


- N 



is? n 2 
n 2 m 2 
-mn mn 


in n mn 


-mn -mn 


J -n 2 


If C ( j is a symmetric tensor, 


-mn mn -n 2 m 2 


C i I * C J i 




22 


12 


21 


( 88 ) 


(89) 


(90) 


(91) 


(92) 


19 
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therefore 


Equation 91 can be simplified as follows: 

C 


C l 2 * C 2I 


'll 


22 


')2 


2 2 

m n 2mn 


2 2 

n m ,-2mn 


■mn mn m^-n^ 


'l l 


22 


(93) 


These are the transformation equations ofT + for second-rank tensors which include the stress 
and inertial tensors. /.• ' 

In + 3-space, the range for the indices will be 1, 2, and 3. If the i-otation remains the same 
as T + before, the following results will be obtained: 


1) When i = j = 1 


C ll ‘ , lm T ln C m« 


, ll , ln C ln + , l2*ln C 2n + f l3 f ln C 3r» 


* * I l (, l l C M + f I2 C I2 + ^I 


+ t l2 t, ll C 2l + t l2 C 22 + 


f ’l3 (, | l°3l + *I2 C 32 1*13*33* 


Since for the present T-matrix, 


(94) 


*13 * *23 * *31 * *32 * 0 


(95) 


then 


® " 2 2 
Cj, * mC ( J + n C 22 + 2m»C |2 


(96) 


where C { j is assumed to be a symmetric tensor, and the results of Equation 96 are the same 
as that in Equation 93. 


20 
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2) When i = 2, j = 2, 


22 


1. 1. C 
2m 2n mn 


C. + C- + t-_t_ c, 
21 2n In 22 2n 2n 23 2n 3n 


2 2 

n C ( j + m Cgg ~ 2mnC ( 2 


This is the same as Equation 93 because of the relations in Equation 95. 
3) When i = 3, j = 3, 

Cl- - t, »_ C 
33 3m 3n mn 


(97) 


*3l*3n C ln + 1 32 , 3n C 2n + *33*3r. C 3n 


(98) 


^33 C 33 * C 33 


where Equation 95 is used and t^ * 1. 

4) When i 55 2, J ■ 3, 

C 23 * ^2m*3«S»s 


5) When i * 3, J * 1, 


* *21 *3s C ls + *22*3n C 2s + *23*3n C 3n 


*2(*33 C I3 + *22*33 C 23 


* - «C |5 4-mC l5 


C 3» * *3m*ln C mn 


* *3l*l» C l» * *32*ln °2n * *33*ln °3n 


*33 ^*11 C 31 + *12 C 32* 


mC 3l + " C 32 


(99) 


( 100 ) 


21 


*4* ■ <Wl| g>4i 


-->VV-r' 

\ W v - 
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6) Wher .i = 1, j = 2. 

C I2 1 , im , 2n C mn 


», ,1, C, 4 t, C, 4 t t_ C, 

It 2n In 12 2n 2n 13 2n 3n 


'.."zi 0 ., ♦ , 22 C , 2 ’* , I2‘V 2 , *’2 2 C 22 ’ 


- -mnC (( 4r,-!: 22 4(m 2 -n Z lC |2 

This is the same as Equation 93. Thus, the transformation equations for 3-space, which are 
equivalent to Equation 93 in 2-spnce, are 



Note that Cg 3 component is invariant, i.e., cj 3 * C 33 , and C g2 and C 31 interact with each 



22 
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C^ 3 and C 31 would be the transverse shear stresses of a plate with Its plane parallel to the 

1-2 plane. The states of stress before and after a transformation of 8 = 90° are shown in 
Figure 3. 

Note that when 8 = 90°, m = 0, n = 1. These results, when substituted into Equation 102 
yield 

C ii * C 22 

C 22 C l I 
C 33 * C 33 

c‘ = -c ,,03) 

... 23 31 

C‘ s c 
31 z3 

C' s-C 
12 12 

The relations in Equation 103 agree with the results shown in Figure 3. The shear components 
shown in this figure represent positive quantities (the shear diagonals lie between the positive 
directions of the coordinate axes). 

CONTRACTED NOTATION 

A further simplification of the indicial notation is possible with the contracted notation. In 
dealing with fourth-rank tensors, the contracted notation reduces the number of free indices 
from 4 to 2 but expands the range from 3 to 9. The number of components, according to 

4 o 

Equation 1, remains at 3 ■ 81 and 9 - 81 for the normal and contracted Indicial notation, 
respectively. But if symmetry properties are introduced, the contracted notation can be used 
to obtain a considerable amount of simplification. 

The fourth-rank tensor of interest now is the elastic stiffness or compliance matrices, 
S ijkl or c ijkl' T1>ey *PP®* r 111 N*® (Boe rail sad Hooke’s laws as follows: 

*11 “ c IJM fl M (l041 

i } - 

•ij * S iJM*M 1,05) 


where <r^ ■ stress tensor, e^ * strain tensor, and both tensors have the following 9 com¬ 
ponents in 3-space: 



•ll 

*tt 

•tl 


*13 

*tl 

•tl 

*31 

*11 

*lt 

•it 


*3t 


• 11 

<rtt 

•tl 






• V, 




S' 

'! 

i 

T 


i 
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The contracted notation is established by arbitrarily replacing the double- index system for 
a second-rank tensor with a single-index system. Each component in Equation 106 can be 
represented by a single-index system from 1 to 9; e.g., as shown in Table II. 

TABLE II 

CONVERSION BETWEEN NORMAL AND CONTRACTED NOTATIONS O 


Normal 

Notation 

Contracted 

Notation 

*1 1 

•h 

*1 

*1 

(Vi 

*22 

*2 

*2 

*33 

*33 

*3 

*3 

*23 

2 *23 

*4 

*4 

*31 

2 *3, 

*5 

*5 

*12 

2 *I2 

*6 

*6 

*32 

2 *32 

*7 

*7 

*13 

2 *I3 

*8 

*8 

*21 

2 *2I 

*9 

*9 ■ C 


In contracted notation, engineering strain is used instead of tensorial strain and Equations 
104 and 105 can be written as: % 

e, ■ SjjCj (107) 

»■, * C|j 0j (108) 

• •. f C 

where 1, J » 1, 2.and 9. In this notation, the range and summation conventions are 

retained. But some modifications in the Interpretation ot the indices must be made. First, 
the range of free indices no longer correspond* to the number of dimensions in space. Sec¬ 
ondly, the tensorial rank no longer correspond* to the number of free indices. Finally, the 
contracted notation is an artificial notation which happens to provide significant simplification 
in the representation of the generalised Hooke’s laws but only second and fourth rank tensors 
where 3-space is involved. The use of contracted notation for other tensorial relations should 
be limited to special cases. ' 

Returning to Equations 107 and 108, the range of the indices can be reduced from 9 to 6 if 
the stress and strain tensors are symmetric l.e., 

<r,j * Ojj, and «|j « *jj (109) 



* 
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which, in contracted notation, means 

<r 7 =Cr 4‘ a 8 = £ V a 9 = £r 6 ( " 0) 

*7 = *4 ’ 6 8 = 8 5 ' e 9 =# 6 (,,,) 

Thus, the symmetries shown in Equation 109 reduce the number of components of and 

f.i im 81 to 36 in 3-space, and from 16 to 9 in 2-space, Equation 6 in contracted notation be¬ 
comes: 

°ll °l 2 ° 16 °I9 

> -v 

Ojj = ° 2 ‘ ° 22 ° 26 ° 29 (. 12 ) 

°6I °62 °66 °69 

.“SI °92 °96 °99_ 

where the replacement of the 4-index system to a 2-index system follows the relationship in 
Table II. If the is a compliance matrix in the generalized Hooke’s law and the stress and 

strain tensors are symmetric. Equations 110 and 111 can be used to simplify Equation 112 
as follows: ; 


:*.. 

*12 

S I6 

*21 

S 22 

S 26 

.* 6 . 

*62 

*66 


The S.. in 3-space will be: 


* 1 . 

* I2 # 

*13 

*14 

*15 

*16 

V 

*22 

, *23 

*24 

*25 

*26 

* 3 . 

*32 

*33 

*34 

*35 

*36 

* 4 . 

*42 

*43 

*44 

*45 

*46 

*31 

*62 

*63 

*64 

*65 

*56 

8 S . 

*62 

*63 

*64 

*66 

*66 


An additional symmetry property can be established from the strain energy consideration. 
If the existence of an elastic potential is assumed then an increment work per unit volume is, 


dW * o-.dSj 


♦V 


25 
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then substitute 


By integrating, one gets 


a i * C 1J»J 
dW * CjjSjdSj 


W - gCy.jS, 


( 116 ) 

( 117 ) 


( 118 ) 


Similarly, one can show that 


W 


■ K'l'i 


(119) 


From the elastic potentials, one can derive the generalised Hooke’s law, as follows: 

dW 


i *Vi 


( 120 ) 


and 




-» Cn 

*•»*•! 1 


(121) 


Similarly, 


fc.fe, 


(122) 


Since the order of differentiation Is Immaterial, then 

0 • • • C11 

•I 1» 


( 123 ) 


i.e., the stiffness matrix must be symmetric with respect to the indices in the contracted 
notation. In a similar manner, one can show that 


S «S 

»J 11 


( 124 ) 


This additional symmetry will simplify Equations 113 and 114, such that there are only 6 and 
21 independent components in 2-space and 3-space, respectively. 


The conversion between the normal and contracted notations for S., and cannot be 

derived directly from Table II. Some correction factors must be introduced because of the 
difference between the engineering and tensorial strains. The factors can be established by 

26 
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. ~ 'r * , \ 

expanding Equations 104, 105, 107, and 108, and by comparing them term by term. For 
example, for i = * = 1, from Equation 104, 

a = S cr + S <r + s O’ 

II Mil l I 1122 22 1133 33 


+ (S II23 + S il32 > °23 + tS t 131 + S l 11 3 )<r 3l 
+ <S I 112 + S l 121 )cr | 2 


(125) 


From Equation 107 


e = S tr + S <r + S o’+S or+S o - + S cr 
I III 122 >33 14 4 15 5 16 6 


(126) 


Comparison of Equations 125 and 126, and an assumption that the matrix is symmetrical 
yields 

S llll * S| l 
S l 122 * S I2 
S l 133 * S I3 
2S .«23 ‘ S I4 
2S II3I * S I5 
2S Mt2 >S l6 


(127) 


By repeating the process similar to that shown la Equations 125 and 126, one can establish 
the following conversions of the components of the compliance matrix:. 


S|JW ■ V f#r • * 1 .2,or 3 

2 S ijkl « Sqr for q • 1,2 , or 3; r ■ 4,5, or 6 
or q>4,5,or6, r» 1,2,or3 


(128) 


• ®or for q,r ■ 4,5,or 6 


qr 


The conversion factors for the stiffness matrix can be similarly established. In 1 * J - 1, 


a M * C ! 111*11 + C 1122*22 4 ' C M3S*33 


+ (C, |88 4 C„ 52 )e 28 4- (C, |S| 4 C, 113)03, 


(129) 


4 (C m2 + C, (2 ,Js |2 


From Equation 56 


<r , * C. .s. + C..i, + C.,1, + C,.*. 4 C,.». ) C, 

I III 12 2 13 3 14 4 15 5 16 6 


(130) 


27 
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From Table II, - 

2 *23 = *4 ’ 2# 3I = *5 ’ 2# I2 = *6 <,3,> 


From the last three numbered equations, 
C 
C 


nil * 

C 11 * 

C 1122 = C I2 ’ 

C 113 3 

- C I3 ’ 

1123 = 

C 14 ’ 

C 11 31 = C 16 * 

C ll 12 

* C 16 ' 


(132) 


By repeating the process for the other relations in the generalized Hooke’s law, one can 
establish that 


C ijkl = C qr 


(133) 


Thus, the conversion factor is unity for all components of the stiffness matrix. But for the 
compliance matrix, relations in Equations 128 must be followed. The contracted notation must 
be handled with care. The tranformation equations must be derived using the normal indicial 
notation. The relation between the components of the compliance and stiffness matrix of the 
two notations must include the proper conversion factors as shown in Equations 128 and 133. 


The use of the contracted notation has often been inconsistent in many current publications 
and reports. In many cases, Censorial strain is retained. In other cases, both single and 
double index systems are used simultaneously, e.g., e^.e^, and e 12 for the strain components 

in 2-space. Finally, e J2 is sometimes represented by e 3 instead of e g . 

The conversion between the normal and contracted notations as listed in Table II and 
Equations 110 and 111 are recommended because a consistent notation between 3-space and 
2-space, and consistent operations (range and summation conventions) can be retained. In 
the contracted notation, fourth-rank tensors are represented by double-index quantities 
(C . S^); second-rank tensors, by single-index quantities. The range for 3-space is 1, 2, 3, 

4, 5, 6, and 2-space range is 1, 2, 6. The latter range is probably better than 1, 2, 3 because 
it avoids the similarity between the range of 3-space of the normal notation and the range of 
2- space of contracted notations. 

h 

Finally, the contracted notation as listed in Table II may be considered authoritative 
because it follows the notation listed in many textbooks (References 1, 2, and 3). 

FOURTH-RANK TENSORS 

The transformation equations for fourth-rank tensors contain 81 equations. With the aid of 
symmetry properties, e.g., ojj * ojj.e^ -e^, and ■ S^ ( and (the last two equations 

are in contracted notation), there are only 21 equations. A further simplification can be intro¬ 
duced if the T-matrix is limited to a rotation about the 3-axis, as shown in Figure 1, whe j, 




(134) 
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The transformation equations must be derived from the normal indici&l notation, not the 
contracted notation, because in the latter notation artificial simplifications are introduced. 
The number of free indices no longer corresponds to the tensorial rank. Repeating the form 
of Equation 30, in which 

C* - ♦ • ♦ | r — 

ijki im jo kelp mnop 


tne transformation of S.^j is, similarly 


S',,., = t. t. t u t, S 
ijfcl im jn ko Ip mnop 


(135) 


Further, in 3-space, by taking advantage of l^ 3 = t r g = t 31 = t^ = 0: 

S' 


11II 


t, t. ♦, S 
Im In lo Ip mnop 


Summing m gives 


*1111 ” *11* lo*lo*l p S loop * *l2*lo*lo*lp S 2oop 

+ *l^ln*lo* lp S 3nop 


Summing n gives 


S T III *11 ( *l l*lo*lp S llop + f l2*lo*lp S iaop + \*{lo*lp S l3op 1 


+ t. s_, t. s„ + 

12 II lo Ip 2lop 12 lo Ip 22op 


Nflo* lp S 23op > 


Summing o gives 

• S' , ■ t 


llll' ’ 11 [ *1 l* t l.l , l.p l 


8 l llp * * l2*»p S M2p +N S^lp i 


8 I ISp* 


* , ia <t U t lp S l2lp + *l2*lp 8 l22* *' > S^lp 8 l23p^] 
+ *12 [*ll (, ll*lp 8 2llp * *!2*lp 8 2l2p *^S^Ip 8 2l5p^ 

+ *l2 <f | I t lp 8 22lp + t l2 , lp*222p 4 ■ > ^»Ijp 8 ^25p , ] 
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Summing p gives 


{ f t 1 [ t ll <t M S lll I 41 ■* > I2 S I 112 + ^>S?!i 13* 


+ s..«, + ♦,-S.i - + > K»S 11 „)1 

12 It 1121 12 1122 r3. 1123 J 

+ f l2 [ f l S 4t l I S I2II 4, I2 S I2I2 4 \^ S I2I3 ) 

+ t l2 (t M S l22t 4 t |2 S l222 +,SV ^ S (223 > ] } 

4 f l2 Oil [ f l I ( *1 l S 2l I I +t l2 S 2ll2 + '*Hf2ll3 ) 

+ t j 2 ( , , | s 2i2, +t (2 S 2l22 4 \^ S 2I23 ) ] 


+ *12 [*11 | S 22I I 4 *I2 S 22I2 4> V^22I3 1 


4 t i2 (t l l S 222l 4 *I2 S 2222 4> >*?2223 , ] } 


4 3 2 2 

nr% S| | || T in nS| | ^ t (n nS| |^| iw ft Sj 122 

, 3 c . 2 ± 2 2 e . 3 e 

4 m nS |2| j + m " *,2|2 4 * n S |22| 4 m " S |222 


4 "Hill 4 




♦ ib n 


+ IBB 


4r " 2 "% 2ll 4 mn %2l2 4 mn ^®222l 4 ,,4s 2222 

Contracted notation, with the proper conversion factors, modifies Equation 136 to 


s; . - IS * nt^/2 ♦ m sl^/2 4 m n S, 2 

+- ib S b 8 #| /2 4 m 2 n 2 ^ 6 /4 4 si 2 i» 2 ^ € /4 4 mn^y^ 
4 mS »*g| /2 + m*sS^ f /4 4 ib^ 2 ^/# 4 


o J t l 

4 m iiB j( 4 mu V 24 mn 


V 2 


From Equation 124, where 8^ is symmetric, Equation 137 can be further reduced to 

S',, * m 4 S,, 4 2m 2 nS, 2 4 2m 8 n8 |S 4 n 4 S 22 42mB S S 26 4 m 2 n 2 S 66 (138) 

By following the steps described above, the transformation equations for the temalning 20 
components of the compliance matrix and all the components at the stiffness matrix can be 
derived. The transformation, like Equation 138, applies to a transformation consisting of a 
proper rotation about the 3-axis. This is a very special transformation. A general transfor¬ 
mation, for which all components of the T-matrix are nonzero will result in transformation 
equations considerably more complicated than Equation 138, in which only 6 out of a total of 
21 components affect the Sj component. The transformation equations for S and C^, when 
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subjected to a proper rotation about the 3-axis, can best be presented in the following tabular 
forms. The use of the tables can be seen by comparing Equation 138 with the first row of the 
first table. In case of the stiffness matrix C^, appropriate factors shown in the column and 

row headings must be properly incorporated, as follows: 


c ii = m4c ii 


2 2 3 

+ 2m n + 4m nC|g 


n% 2 + 4mn 3 C 26 


2.2 

+ 4m nC„ 


66 


(139) 



W 

W 

S ,6 (2 V 

S 22 (C 22 J 

S 26 (2 S 6 ) 

S 66 (4C 6 6 ) 

S V, (C V.> 

4 

_ 2 2 

„ 3 

4 

■' 3 

2 2 

m 

2m n 

2m n 

n 

2mn 

m n 

S I2^ C I2^ 

2 2 

4 4 

3 

2 2 

3 

2 2 

m n 

m + n 

mn_ 

-m n 

m n 

m n_ 
-mn s 

-m n 

s ; 6 (2c, «6 ) 

_ 3 

_ 3 

4 

„ 3 

- 2 2 

3 

-2m n 

2m n 3 
-2mn a 

m 2 2 
-3mn 

2mn 

3m^n 
- n 

m n_ 

-mn 3 

S 22 <C 22 } 

4 

- 2 2 

« 3 

4 

. 3 

2 2 

n 

2m n 

-2mn 

m 

-2m n 

m n 

S 28 (2C i«> 

-2mn 3 

2mn® 
-2m n 

3m 2 n 2 

-n 

2m 3 n 

i . 

4 2 2 
-3m n 

mn? 

-m n 


2 2 

2 2 

3 

2 2 

3 

. 2 

4m n 

-8m n 

4mn_ 

-4m i* 

4m n 

4m n_ 
-4mn 

(m 
-n ) 


(140) 



S »S <<5 IS» 

w 

S 36 ( *W 

s ;v c ;3» 

[■ : y m 2 , 

.2 

II • 

mn 


n 2 - 

. 

• 2 
m 

-mn 


-2mn 

2 mn 

2 2 
m - n 


V C 44» 

W 

w 

S i4 (C M» 

m 2 

-2mn 

n 2 


mn 

2 2 
m - n 

-mn 

S « (C i 8 ) 

2 

n 

2 mn 

2 

m 


(141) 


(142) 
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rj • 

VV 

W 


m 

-n 

S 3^ C 35 > 

n 

m 


(143) 


S ’«4 { SV 
s, . 5 tc ; 5 > 
S 24 (C iV 
^ 5 ( C 25 ) 
S 46 l 2C 46 ) 




S !4< C «4> 

V C I5> 

S WSt> 


W 2 <W 


3 

2 

2 

3 

2 

2 

m 

-m n 

mn 

-n 

m n 

-mn 

2 

3 

3 

2 

2 

2 

m n 

m 

n 

mn 

mn 

m n 

2 

3 

3 

2 

2 

2 

mn 

-n 

m 

- m n 

-m n 

mn 

3 

2 

2 

3 

2 

2 

n 

mn 

m a 

m 

-mn 

-m n 

2 

2 

2 

_ 2 

3 

3 

-2m n 

2mn 

2m n 

-2mn 

m 2 
-mn 

n 2 
-m n 

- 2 
-2mn 

-2m 2 n 

2 mn 2 

2m 2 n 

2 

m n 
-n* 

3 

m 2 
-mn 


(144) 




(145) 


Note that under a particular tranafoimation, which In this case la a rotation about the 3-axis, 
the components of and art) arranged in 6 groups, each of which is numbered as an 

equation betw-en Equations 140 and 146. The components of each group will Interact with one 
another, but a *e completely uncoupled from the other groups. This information is useful in 
the study of t laatlc symmetry. If for a given material all the components within a group are 
zero, they wtd remain sera for all angles of rotation about the 3-axis. 

'T ' V'.. 1 • 

INVARIANTS ’ ■ 


There are a number of invariants associated; with the compliance and stiffness matrices 
with respect to the rotation about the 3-axis: 


I, ♦ S-- 4 28,- 

I VII 22 12 


From Equation 140 


- , 4 . 4 , _ 2 2,_ 22._ 2 2 4„ 4,- 

Ij » (m 4 n 4 2m n )S t ^i ; 4- <2m n + 2m n 42m 42n )S j2 

+ (2m 3 n -2mn 3 4 2mn 3 2m 3 n)S j6 4- (n 4 4 m 4 4 2m 2 n 2 )8 22 

+ (2mn 3 - 2 n ?n 4 2m 3 n - 2mn 3 )S 26 4 (m 2 n 2 4 m 2 n 2 - 2m 2 n 2 )S, 


66 


(146) 


S 11 4 S 22 4 2S I2 
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*2- S 66- 4S .2 

9 9 ' 7 9 44 22 

(4 m n- 4m it-lS - 4(m 4-n 4 2m n )S (2 

3 3 33 22 22 

4 4(mn - m n - mn + m n ®|g + (4m n ~ 4m n )S 22 

3 3 ° 3 : 3 T 2 2 2 2 21 

4 4(m n - mn - m n 4 mn )S^g +• [ (m - n ) 4 4m n JS, 


66 


(147) 


, \. s 

• 4 


S 66 - 4S .2 


h = S 44 * S 55 

= (m 2 4 n 2 )S 44 4 (2mn - 2mn)S 45 4 (m 2 4 n 2 )S 55 


(148) 


S 44 +S 55 


V S 23 +S I3 

S (m 2 4 n 2 )S )3 4 (m 2 4 n 2 )S 2J 4 (mn - mn)S 36 


(149) 


S .3 +S 23 


V S 34 + ^V 


2 2 2 v 2 2 2 

(n» 4 n 4 (r> 4 n ■Sb 


s 2 ..2 

S 34V S 38 


t-- 

V; 

■1 


2 « * S 33* *33 


(150) 


(151) 


An anisotropic body is subjected to hydrostatic pressure so that 


and 


“ * 9 *' *% * P 


°4 “ a 8 s U " 0 


The change in volume is the sum of the normal strains; i.e.. 


8 v / v = s , 4S 2 4. 3 


I 

f 
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The shear strains e^, e^, and eg under hydrostatic pressure will not be zero for an aniBo- 

troptc body, but they do not contribute to any change in volume. Adding the first three equations 
of the generalized Hooke's law, 


;; .. *i = 

'irt + V!*Vj 

i *2 = 

S l2*i ♦ 

I ;i. * 3 = 

S |3 tT | + S 23 a 2 * S 33°3 

8v/v = 

[ S M * S 22 + S !3 + 2IS (2 


Thus, the compressibility of an anisotropic material is 

,v/8v 1 / [s„+s 2J ts„+j(s |2 + s 23< s 31 i] 


1/(1, +I,«I 4 ) 


(152) 

(153) 


The compressibility Is also invariant. Similar invariants for the stiffness matrix can be 
established immediately, as follows: 


C M + 

C 22 * 

2C* 

12 

• 

C .. + 

C 6« ~ 

C l2- 

c e« 

• 

C I2 

C i4 + 

C 3B " 

C 44 

+ 

C 55 

C » + 

'IS¬ 

C 23 

+ 

C t3 

C i4 2 + 

C' 2 

C 35 

* C 34 2 

+ C 38 ! 


“I 

J 2 

J 3 

J 4 

J 5 


J « * C 33 “ C 33 


'22 


'12 


( 154 ) 
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SECTION ID 

ELASTIC SYMMETRIES AND ENGINEERING CONSTANTS 


The compliance and stiffness matrices in 3-space contain 21 independent components. If 
a material contains symmetry properties with respect to certain directions, which can be de¬ 
scribed in terms of coordinate transformations, the number of independent Components will 
reduce. Ultimately, if a material is isotropic, there are only two independent constants. In 
this section, a few commonly encountered material symmetries and the relations between 
S.j or and the engineering constants for various materials will be examined. 

- IN¬ 

ELASTIC SYMMETRIES 


A triclinic material is the most general anisotropic material where all 21 elastic constants 
are independent. A fourth-rank tensor in 3-space will have 81 components. If both stress and 
strain are symmetric tensors, the 81 components can be represented by 36 independent com¬ 
ponents. If the stiffness and compliance matrices are symmetric, which can be demonstrated 
by assuming the existence af appropriate elastic potentials, then only 21 of the 36 components 
are independent. 


Further reduction in the number of Independent components can be introduced if additional 
material symmetry exists. One of the simplest forms of symmetry is the monoclinic material 
which possesses one plane of symmetry. Let plane x-y, or equivalently, z - 0, be a plane of 
symmetry; then the properties at +z are equal to those at -z. If a coordinate transformation 
of 

r« o oi 



I » 0 


(155) 


Lo o -i J 


will leave the and Intact, the material by definition is a monoclinic material. The stress 

and strain components will transform .with results very similar to those shown in Equations 
94 through 101. Both the stress and strain components after the transformation will be as 
follows: 

C i l “ C ll • C 22 * C 22 ' C S3 * C 33 ' 



-C 


23 ' 



■ -C 


31 ’ 




where C^, following Equations 94 through 101, is a typical second-rank tensor, and is not the 

stiffness matrix in the contracted notation. In the contracted notation, only the following com¬ 
ponents of stress and strain change signs: 


°4 * ~ **4 ' 

V s *4 *5 * *5 


(156) 
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Compare the first equation of the generalized Hooke’s law written in the new and old 
coordinates: 

"i = C . + C I2*2 + C .3«3 + C .4% + C | 5 «5 + C ,6*6 


cr { ~ C l| « ( 4 c , 2 « 2 +c , 3 » 3 ~ c l4 ® 4 -C, 5 e 5 + C |6 e 6 


If cTj = <Ty it is necessary that 


C .4 =C I5 = ° 


(157) 

By considering the remaining five equations in the generalized Hooke’s law, one can show that" 


C 24 = C 25 “ C 34 “ C 35 = C 64 = C 65 * °‘ 


(158) 


The compliance matrix must allow the same components to vanish as those shown in Equations 
157 and 158. The number of independent components for S. and C.. reduce from 21 to 13 for 
monoclinic materials, as follows: lJ 




C 



S I2 ®l 3 

0 

0 

S .6 

S 22 

0 

0 

S 26 

S 33 

0 

0 

*36 


S 44 

S 45 

0 


*65 0 


C 22 

C I3 

0 

0 

C I6 

C 22 

C 23 

0 

0 

C 26 


C 33 

0 

0 

C 36 



C 44 

C 46 

0 


(159) 


(160) 


66 


Since both S.. and are symmetric, the lower half of the matrices is not shown. The com¬ 
ponents shown above correspond to a monoclinic material with z = 0 as the plane of symmetry. 
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If a different plane of symmetry exists, say x = 0, the nonzero components of S.. and C.. will 
be different from those shown above. The matrix will be: ■ * * 


S |2 S l3 

S I4 

0 

0 


S 22 S 23 

S 24 

0 

0 


S 33 

Ha 

0 

0 

(161) 


S 44 

0 

0 




S 55 

S 56 





S 66_ 



The number of independent components for all monoclinic materials remain at 13 irrespective 
of the orientation cf the symmetry plane. 


If an anisotropic material possesses two orthogonal planes of symmetry, say, x = 0 and 
z = 0, the independent components must satisfy the in Equations 159 and 161 simultaneously. 

This material is called orthotropic and must have the following compliance matrix: 




S I2 

S I3 

0 

0 

0 

^2 

S 23 

0 

0 

0 


S 33 

0 

0 

0 




0 

0 


\- 

44 

Sb 

0 


86 _1 



(162) 


The number of independent components reduces from 13 to 9. If a material has two orthogonal 
planes of symmetry, it will automatically have symmetry with respect to the third orthogonal 
plane. 


If a material has a plane in which the property is Isotropic, this is called a transversely 
isotropic material. Let us assume that the x-y (or 1-2) plane is isotropic, there is no pre¬ 
ferred orientation in this plane. All properties will remain invariant under a rotation about 
the z-axis. The indices 1 and 2 in the and are Interchangeable, thus: 


II 


* S. 


22 


II 


• C, 


22 


S « S 
13 23 


: I3 * C 23 
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'II 


'12 


13 


'13 


'33 


0 

0 

0 

C 44 


0 0 

O 0 

0 0 

0 0 


44 


( C ,|- c , 2 )/2 


(166) 


If the isotropic plane is in the 2-3 plane or with respect to the 1-axis, the compliance matrix 
becomes 


’ll 


S I2 

S ,2 0 

0 

0 

^22 

S 23 0 

0 

0 


*22 0 

0 

0 


2( V- S 23 ) 

0 

0 




S, 


'55 J 


(167) 


The number of Independent elastic constants for this material is five. 

In the case of isotropic materials, indices 1, 2, and 3, and 4, 5, and 6 are interchangeable; 
thus. 


hs m hi m V’ 

*55 ‘ *•« • 


8^- 


11 


C = C 

22 33 


C »* C 3I* C I* 
C «4‘ C »» * C M 


In addition, 


8 44* 2(, l I ’ *12 1 ’ 


C s (C • C ) /z 
^44 "'ll W I2' 
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The final S.^ and for isotropic materials are 


J I1 


S t2 S I2 


II 12 


II 


0 

0 


0 

0 

0 


2(S — S I 0 
II 12' 


0 

0 

0 

0 


2(S - S ) 0 

II 12 


2tS ..- S .2 > J 


( 168 ) 


V 


'll 


'12 12 


II 12 


0 

0 

0 


0 

0 

0 


(C,,-C, 2 >/2 0 


0 

o 

o 

o 


{C .,- C .2 ,/2 0 


( V 


C .2> /2 J 


(169) 


There are only two independent elastic constants. 

ENGINEERING CONSTANTS 

Engineering constants usually refer to Young's moduli, Poisson's ratios, and shear moduli. 
These material constants can be measured from simple tests, such as the uniaxial tension 
or pure torsion tests. The constants are therefore more familiar and understandable than 
the components of S^ and C^. The relationship between these components and the engineering 

constants will be established in this section. 

Since most simple tests are performed with a known imposed load or stress, the resulting 
displacement or strain is measured. The former is the Independent variable; the latter, the 
dependent variable. Thus, the components of the compliance matrix can be more explicitly 

determined than those of C^. The following relations between the components of and the 
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. 

c* 

S 

' 

engineering constants can be established Immediately from (he nature of uniaxial and pure 
shear tests: - . ' 

S l 1 ,/E »l ’ S 22 = l/E 22’ S 33 =I/£ 33 

- ■- < 

5 

S .2 - 

~ V t2 /E l 1 ’ 

S 23 * '"a' E 22 ' 

S 3I = ■ ,/ 3l /E 33 


1 

= 

~ V 2\ /E 22 


” V | 3 7 E l 1 

(!70) 


S 66 = 

1/6 .2’ 

S 5s‘ l/S 3 ' 

S 44 S,/G 23 


i 

S I6 

1, I6 /E I 1 ’ 

S 26 5 1, 26 /E 22’ 

S 36 = 1, 36 /E 33 


i 


I 


I'.,. 


where the rj -j are the shear coupling ratios. Other components of S,^ can also be expressed 

in terms of engineering constants. But new and unfamiliar engineering constants must be 
invented and it is doubtful that a useful purpose is served by forcing a complete equivalence 
between the engineering constants and the components of S.j. In fact, only orthotropic, trans¬ 
versely isotropic, and isotropic materials are being investigated in these notes. The triclinic 
and monoclinic materials are being omitted to avoid using unfamiliar engineering constants. 
The shear coupling ratios often appear in two-dimensional problems. 

For an orthotiopic material, one can express the components of in terms of the engi¬ 
neering constants. This can be done in a straightforward manner* by making a substitution of 
Equation 170 into un inverted S^. The resulting relations are listed below: 

C l. " (l -*'23*W VE «. 

C 22* ~ r 5l *13**28 


C I2 " <r 2l + *23 r 3l ,VE lI * *12 +>, 'l3 , i2 ,VE 22 
C 13 * ^31 4 *'2l l '32 ,VE n " *I3 + *'23 , '|2 ,VE 33 
C 23 * <*32 +V i2 K 3l ,VE 22 8 <*23 4 *'|3 ,VE 33 

C 44 ‘ °23 
C 3B*°3I 
C «« * ®I2 


( 171 ) 


where 


V : ( I • V V -V V 

v ■ 12 21 23 32 


' V 3l V t3 ' 


( 172 ) 


The relations between and the engineering constants are considerably more complicated 
than those for S^. Poisson's ratios are responsible for the complicated relations above. If 


i. . 

I 
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all Poisson’s ratios are rem which means that there Is no ooupling between the normal strains 
or stresses, then 


C M “ E «.» 


C 22 * E 22 ’ C 33* E 33 1 


C I2 SC I3 =C 32 = 0 - 

If the material is transversely isotropic with plane 2-3 as the isotropic plane, the com¬ 
ponents of Equation 167 can be related to the engineering constants as follows: 


S|, = l/E |{ , 


S 22 = S 33 S ' /E 22 


S I3 = S I2/ “ V |2 /E l I = _I 2| /E 22 

S 23 S “ y 23 /E 22 

S 44- 2(,+ V /E 22 


The components of are 


»VE„ 

C 22 “ Ss* (l " l 'l2*2l ,VE 22 

C I2 * C .3^2. n+ W VE ..*V ,f V VE 22 

C 23 , %3 + ".2 , '2. We 22 


C 44- (C ! 


22 


C 2J WZ * < 1 ‘*23 ’ * V \Z\\ m iz' 2 


C 38 * C ««- °I2 


(173) 


(174) 


whore 


^ 9 I 

V a ( I - 9u y — V - 2V V V ) 

1 *12 21 23 12*21 2S 1 

-I 


[ (I (< “ *’23 


«*] 


(175) 


Simplification of Equation 174 can be made by using Equation 175. If the material is isotropic, 

S,, = S = S = I /E 
II 22 33 

(176) 


S I2 =S I3 = S 23 S ‘* /E 
S 44 =S 55 =S 66 S 1/6 s2(, + *' ,/E 
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The components of C.^ are: 


(I -v)E 

C 11 = C 22 = C 33 " (l + v) (I - Zv) 

V*- 

C |2 = C 13 S C 23 " (1 + v)(l-2iO 


(!77) 


C = C - C = G = --- 

4 4 5 5 66 2(1 + v ) 


In this section, the components of and C.^ for orthotropic, transversely isotropic, and 

isotropic materials are expressed in terms of commonly encountered engineering constants. 
The components of S^. have simpler relations than those) of C.. with the engineering constants. 


TRANSFORMED S tj 

The nonzero components of for monoclinic, orthotropic, and transversely isotropic ma¬ 
terials in coordinate systems other than their material symmetry axes can easily be estab¬ 
lished from the transformation equation listed in Equations 140 through 145. The transformation 
being investigated is restricted to a proper rotation about the J axis. All components which 
have primes designate transformed components. 

In the monoclinic case with the symmetry plane contain. % the 1-axis and the 2-axis, 
Equation 159 represents the independent components of S... A ro ition about the 3-axis would 
result in the following nonzero components: J 


V, 


II 


12 

*22 


*15 

'25 


'55 


0 

0 

0 

®44 


0 

0 

0 


’16 


S' 


Si 


45 


'55 


26 

S 56 

0 


'66 J 


( 178 ) 


The number of nonzero components (20) of Equation 178 does not differ from that of Equation 
159 because the axis of rotation coincides with the normal to the symmetry plane. The 

in Equation 161, however, has a plane of symmetry containing the 2-axis and 3-axis. When 
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the same transformation as above is applied to this monoclinic material, the resulting S. 
becomes r • 




S »i S 


'II 


S I2 

8 13 

8 14 

*;»■ 

S '.6“ 

S * 2 

S 23 

S 24 


S 26 


S 33 

.*84 

S 35 

8 36 




** 

S 46 




S 55 

S 56 





S 66_ 


(179) 


There are 36 nonzero components, as compared with 20 in the principal direction of S.^. But 

this is not a triclinic material because of the 36 transformed components, only 13 are inde¬ 
pendent. To distinguish between the Sj. in Equations 159 or 178 and 179, the former may be 

called special monoclinic, the latter, general monoclinic. The special monoclinic refers to 
the S.j in its principal axes or the material symmetry axes. All monoclinic materials have 

only 13 independent components; and only the special monoclinic material has 20 nonzero 
components, as shown in Equation 159 or 178. A special orthotropic material is shown in 
Equation 162. Using the transformation Equations 140 through 145, the general orthotropic 
material can be shown to have the following S^: 


(180) 


s • :■ 

12 

8 I3 

O 

0 

S ,6 

S 22 

S 23 

0 

0 

S 26 

} • 

i < j [ 

^3 

0 

0 

S 36 

5 i 
if 

!l 

l 


JA4 

S 43 

O 

?. 

i: 


- ' 

S 55 

O 


. , ' i! "66 J 

There are 20 nonzero components, of which 9are independent. A general orthotropic material 
has the appearance of a special monoclinic material. 

i; ' 

A special transversely isotropic material relative to the 3-axis is shown in Equation 165. 
For this material, i< 

S ll * 8 22 * *66* 2(8 11 " S I2 ] 

8 I3 S 23 1 8 44* *88 

By substituting these conditions into the transformation equations in Equation 140,141, and 
142 one can show that 


W 8 *’ 8 *" 0 


44 
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The general transversely isotropic material corresponding to Equation 165 is 



Since the rotation is in the isotropic plane, Equation 181 is the same as 165. An analogous 
situation occurred between Equations 178 and 159 for the monoclinic material. 


For a specia transversely isotropic material with a different isotropic plane, e.g., the 
2-3 plane as sh wn in Equation 167, the corresponding general transversely Isotropic ma¬ 
terial will be 



(182) 


There are 20 nonzero components that are similar to the special moooclinlc material, but 
the number of independent components for a transversely Isotropic material remains at 5. 

For isotropic material, one can substitute the relation in Equation (176) and readily show 
that the remains same as Equation 168, which has 12 nonzero components (2 independent 

components). / 
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” SECTION IV 

TWO-DIMENSIONAL COMPOSITES 


Two-dimensional formulations are of particular interest in the study of composi te materials. 
In this section, the assumptions of plane strain and plane stress, which represent two special 
two-dimensional problems, will be described. Laminated composites as special cases of 
plane stress will also be discussed. All two-dimensional problems are reduced from a three- 
dimensional special monoclinic material with z = 0 as its symmetry plane. Any rotation about 
the z-axis will not make the special monoclinic into a general monoclinic mrterial; this is 
shown in Equation 178. The use of this material which possesses 20 nonzero (13 independent 
components will include, as special cases, the general orthotropic and general transversely 
isotropic materials, shown in Equations 180 and 182, respectively. 

PLANE STRAIN - 


(183) 


\ state of plane strain is obtained by assuming that 

s «/*,. * 2 > 

“i' V*. • 's' 


U_ « u- 
3 0 


where a are the displacements along 1, 2, and 3 axes and u Q is a constant. From the strain 
displacement equation, one obtains 


When this Is expanded, 



*J.» 




* m 3 


■ 0 


4 S3 fcg tog 



O 

(164) 
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For a monoclinic material with the following compliance and stiffness matrices: 



S l( S 12 

*13 

0 

0 

S “I 
5 I6 


S 22 

^23 

0 

0 

S 26 

s.. = 

V 

S 33 

0 

0 

S 36 

1 J 



S 44 

S 45 

0 





S 55 

0 


- 




%6 _ 


” C M c !2 

C I3 

0 

0 

C 16 ~ 


C 22 

C 23 

0 

0 

C 26 

c iJ ■ 


C 33 

o 

■a o 

■a 

0 

C 45 

C 36 

0 





C 55 

0 






C 66 _ 



(185) 


(186) 


it can be shown by substituting Equation 184 into the generalised Hooke’s law in terms of the 
stiffness matrix that 


‘W 0 

°3 * St*l + C 3Z*2 + C 36*6 


and from the generalised Hooke's law in terms of the compliance matrix, 

+ S 32*2 + S 3»°6 > 

Since <r 3 is now dependent on the other stress components, it can be eliminated from the 
generalized Hooke’s law. The results are: 


where 



Vi 



Vis 

S„ 


i,| * I, 2, 6 


(188) 
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■ ' 

Rjj represents reduced constants derived from S^.Itis applicable for the case of plane strain 
Imposed on a special monoclinic material. The components of remain unchanged, so that 

the generalized Hooke’s law for plane strain is 

\ 

°i=V) U89) 


where t, j * 1, 2, and 6. Where i, j - 3, 4, and 5, the Hooke’s law is expressed by Equations 
187. 


For a special orthotropic material, under a state of plane strain, the compliance and stiff¬ 
ness matrices can be written down from the results shown in equation 162, so that 





= 0 


Thus, 


16 



“ 0 . 




V 



S I3 S 23 


33 




0 


0 


(190) 



66 




(191) 
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For a special transversely isotropic material shown In Equation 167, 



Finally, for isotropic material, 
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The reduced compliance matrix In terms of the engineering constants can be obtained by 
direct substitution of Equation 170 into 190, 192, and 194, and the results are 

1) Orthotropic material: 


R ll * 

11 ” "l 3*'3l* /e i 1 

R 22 = 

11 • 

jo 

M 

n. 

-‘V* Vsz’' 6 )) 

r 

- + Va v hz 

R 66- 



'1961 


2) Transversely isotropic material: 


R lt = (l l/ l2 V 2l ,/E tl 

*«■ "-4 ,/E 22 

R |2 * " ’i** 1 +, 23 ,/6 ll 
* ‘W * *23 I/E 22 

V* W' ,e i 2 


(197) 


3) Isotropic material: 


R ll * **22 “ (l ' V ,/E 

R| 2 * -Ml 4 M/E 

%e" 2(1 


(198) 


The stiffness matrix in terms of the engineering constants are exactly the same as those 
for the three-dimensional case. For the special orthotropic material, they are shown in 
Equations 171; the special transversely isotropic material Equation 174; and the isotropic 
material. Equation 177. 


It can be shown that 


Vl« 1 *i D 


(199) 


which indicates, as expected, that the reduced compliance matrix is the inverse of the stiff¬ 
ness matrix. 
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PLANE STRESS 


Another two-dimensional problem can be formulated by assuming 

<r = <r - <r - 0 

3 4 8 


(200) 


From the of the spec ial monoclinic material shown in Equation 185, 


A : A S 0 

4 5 

- S_,a + S or + S_ cr 
3 311 32 2 36 6 


( 201 ) 


The generalized Hooke’s law becomes 


t. = S..<r , whara i,j = I, 2, 6 


( 202 ) 


From Cjj shown in Equation 186 


I 


*3 = ~C 33 ,C 3I*I + C 32*2 f %6*S ) 


(203) 


The e 3 is not an independent component. This is analogous to oj being a dependent component 

in the case of plane strain. Substitute Equation 203 into the generalized Hooke’s law in terms 
of Cjj$ 'then ,-c " 

(204) 


<Wi ;V It;. : 


where 


• » . ;r\. 
U ■ i 
:■ W- 


A.. • raduead stiffness matrix 
C C ’ - 

i3 13 

* c «i - “cl 


(205) 


’33 


and 


i.J ■ I, 2, • 


Thus, for plane stress, the remains the same as a three-dimensional material, whereas 
the must be replaced by Qjj. 


For a special orthotropic material, 


S «S « S « S >0 
16 26 36 45 

C Id * C 28 * C 36* S 45 s 0 


(206) 
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The and for a state of plane stress are 



For a transversely isotropic material, 



% 2 * 8 55 * 

S I 2 * S 13 


1 

C * C 

22 35 , 

C ■* C 

12 13 


.s 

J 

& 

Then, is the same as Equation 207 , but 

is 


i 

\V 

r c . 

n cTT 

22 

C I 2 C 2 S 

0,2 ' <*2 

0 



Finally, for an isotropic material, 


S H 

CM 

CM 

CO 

it 

S l2 * s 23’ 

S * 

66 

2(S | | “ S I2 1 

c .l 

Si 

o 

•1 

C I2 * C 23 ’ 

C * 
66 

<C .I - C I2 ,/2 



(207) 


(208) 


(209) 


( 210 ) 


( 211 ) 
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Then, Sjj and are 


S ii = 


S|2 


2(S ..- S .2 , J 


( 212 ) 


C n - C 


1L 

ii 


'12 


12 

C .f 


C ll ~ C 


12 


11 


(213) 


<W 2 J 

In terms of the engineering Constanta, the for all three materials follows Equation 170. 

The can be obtained by substituting Equation 172 into Equation 208 for the orthotropic 

material; Equation 174 into 210 for the transversely Isotropic material; and finally. Equation 
177 into 213 for the Isotropic materia.. The resulting relations are listed below; 


1) Orthotropic materials 
Q 


« E /(I -» v ) 
II II 12 21 


°22* E 22 /0 “ r l2 r 2l* 

Q 12 * , kl B i I /C ' -*2 V**i2 E 22 /( ' ' Vt! 1 


(214) 


0 ■ 0 
*66 12 


2) Transversely isotropic materials; 


Q ll * ^l 7 * 1 " ''ll*21 * 


®22* *22^* ~ r \2 v 2\ * 

°I2 'uW-'tz'z,* 


(215) 


Q r 6 

66 12 •' 

The Q.. for the transversely isotropic material is the same as that of the orthotropic material, 
-ngineering constants associated with the 3-direction, e.g., and v^, do not appear in the 


53 




AFML-TR-66-149 
Pt II 

plane stress cases. These constants, however, appear In the plane strain case. Thus, under 
a state of plane stress, orthotropic and transversely isotropic materials are Identical, hut 
under plane strain, and three-dimensional problems In general, these two materials are of 
course different. 


3) Isotropic materials: 


0 I| SQ 22 SE/{, " V ' ’ 
Q = VE/{I - v Z ) 

Qgj » E/211 + v ) 


( 2 ! 6 ) 


Again, it can be shown that 


Q S = 8 
ij Jk ik 


(2l r ) 


which indicates that is the inverse of S^. 

COMPARISON OF PLANE STRAIN AND PLANE STRESS 

In two-dimensional problems, modifications to the compliance and stiffness matrices may 
be necessary. The appropriate Hooke's laws are 


1) Plane Strain: 


"■* Vi 


Wi 


(2161 


where 


• S { j - 


S i3 S j3 

S 53 


2) Plane Stress: 


where 


U • 


9 \ * Vi 


■* ■ Vi 


( 219 ) 


C I3 C J3 

Si-'M-c— 


i.i * », 2, « 
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In microscopic mechanics analysis of composite materials, it is a common practice to 
solve an inclusion problem in plane ctrsiu arl plane stress, from which the effective trans¬ 
verse stiffness E ^. e.s., can bepredicted. Let us assume that a transverse load o^ is im¬ 
posed, while <Tj - cTg = 0. In a plane strain case, from Equation 218 


*2 = 


Thus, the effective stiffness in the transverse direction (the 2-direction) is 

V* * ,/B a 

In terms of engineering constant?, from Equations 196 and 197, 

W'se'"-’; n'sJ 


* E / (I - J 12221 

22 23 

for orthotropic or transversely isotropic materials, respectively. For a plane stress case, 
with <r 2 / 0, CTj = = 0, 

Ol/e * l/$ « E„ (223) 

2 2 22 22 

This is true for both orthotropic and transversely isotropic materials. A comparison of 
Equations 221 or 222 with Equation 223, shows that the Poisson ratio associated with the 3- 
direction enters the plane strain analysis but not the plane stress case. 

LAMINATED COMPOSITES 

laminated composites to be considered in this subsection consist of layers of thin ortho- 
tropic plates bonded together. Eaoh layer may have arbitrary thickness and orientation of its 
material symmetry axes. In general, each constituent layer is a genera! orthotropic material. 
Assuming that the coefficients of thermal expansion or contraction are also orthotropic, the 
three-dimensional generalised Hooke's law for each layer may be modified as follows: 


• ( * S,|«j ♦•jT , i|| « I, 2, ... , 6 


(224) 


where a ( * thermal expansion matrix, T ■temperature. The first term represents mechanical 

strain, and the second, thermal strain. The thermal expansion matrix is a second-rank tensor, 
as indicated by the single index in the contracted notation. For an orthotropic material, the 
independent components of a second-rank tensor are, for example, j, a^, and a 3 , while 

a 4 = a 5 = a g * 0. In a general orthotropic material with a rotation about the 3-axts, a g will not 

be zero. Equation 224 can be inverted to have the following form: 


a . 'VS " 0 j T) ’ }>i V'*• 


(225) 
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If we aseume that each constituent layer Is under a state of plane stress, which Is reasonable 
for thin plates subjected to ln-plane (plane 1-2) stresses, 

<r 3 = a 4 S<r 5 s 0 ( 226) 

From Equations 224 and 225, we can show that 

*4 = * 5 ' 0 


•»'* T *■ 4“' 


C C 

•, T > • -V'-C^W’ 12271 


From this, reduced constants can be obtained, following a similar derivation described in 

the last subsection. The generalized Hooke ; s law in a state of plane stress including the 
thermal effect is: 


= ® i " OjT) , i, j 5 1, 2, 6 


(228) 


where 


C IJ C 13 

°'i ■ ■ -v 


(229) 


and their relations to engineering constants are shown in Equation 214 or 215. 


If a lamina .' composite is thin and the deflection of the composite plate is kept small 
relative to it thickness, it Is reasonable to assume that normals to the middle surface are 
nondeformabh . With this assumption, 


Sj » s. + zkj , 


» « I, 2, 6 


(230) 


where e” * In* plane strain and * curvature with the following definitions: 

k. * m. 


s° ■ u. . 

1.1 


3,11 


e° » u 
2 2,2 


k ■ it 

2 3,22 




k * 2u 
« 3,12 


To be consistent with the use of engineering strains, both the in-plane strain and curvature 
are expressed in engineering quantities, as opposed to tensorial quantities. This is not a 
trivial point. The use of engineering strains and curvatures is preferred, particularly in the 
contracted notation. A number of symmetries, e.g., in C^, S^, and Q^, and later in 

Ajj, B.j, and in 'Equation 236, can be maintained with the proper use of the contracted 

notation. As stated earlier in these notes, some authors have employed their own contracted 
notation and may have caused some unnecessary confusion. 
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The substitution of Equation 230 into 228 yields 

crj = + zkj) - OijfljT (231) 

In a thin homogeneous or laminated plate, it is convenient to deal with stress resultants and 
stress couples, which are defined as follows: r ' ■ r 

h/2 

(N, , M j) * f <r, (I, z jdz (232) 

-h/2 

The Integration of Equation 231 gives 


where 


n i + n T* v! + Vi 
“i * m i + »*7 * v° 4o ij“j 


h/2 

/ Q,. a.T (I, z) dz 

-h/2 


h/2 


^VVV “ / « u < 


-h/2 


(233) 

(234) 


(235) 


(236) 


Equations 233 and 234 are the basic constitutive equations for laminated anisotropic plates 
subjected to small deflections. These equations are applicable to thin shells, if the radii of 
curvature of the shells are large in comparison with the shell thickness. The effect of tem- 

'• T T 

perature is taken into account by the equivalent thermal loadings, and M ( . The deforma¬ 
tions Induced by a temperature change are equal to those produced by applying the thermal 
loads. ’ 


The stress at any point in a laminated anisotropic body can be obtained by inverting 
Equations 233 and 234, which in matrix form, are: 



(237) 
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at:] 


where 


A* * A 


-1 


B"* * -a"'b 

-I 


H* = BA 
D^s D - BA -, B 
A* s A* - BFD*" 1 H* 
B‘s B*D*~ 1 


0 * 0* 


l-I 


The substitution of Equation 239 into 230 yields 


S, • «, 4-xh, 


» (A'|j 4 iBjj) Nj 4(b|, 4 xO||)M| 


(239) 


(240) 


(241) 


From Equation 231, the stress components for the k-th layer of a laminated composite are 

U42) 


jrjki , ofj'Nsj - oj k hr) 


> 4( V 4 *°i‘k ■, T J 

- CK**v,J 


(243) 


The last equation is derived by assuming a constant temperature in Equation 235. In this 
case the temperature effect is lumped Into one term. 

s' 

Matrices A, B, and D are the intrinsic properties of a laminated composite. They depend 
on the properties of each constituent layer Q^, the thickness h^, and the stacking sequence 
of the layers. 
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If all layers are quasi-homogeneous, the integrations of Equation 236 can be replaced by 
summations, as follows: 

1 ' 1 •' ‘ 

*„ * - s' ,244 > 

B 'l * 

°i| *T,t° !""?♦ 12461 


Thus, matrices A, B, and D are simple to determine for a laminated composite with a limited 
number of layers. But for the determination of stress and strain, from Equations 242 and 241, 
respectively, the prime matrices A*, B*, and D' are required. These matrices are obtained 
by the matrix inversion operations shown in Equation 240. The inversion of a 6 x 6 matrix is 
very difficult *o do by hand. This unfortunate situation is unavoidable in the case of a general 
laminated composite where all 36 components are nonzero. 


A considerable amount of simplification is possible if the B-matrix is identically zero. 
This occurs if the laminated composite is symmetrical with respect to the middle surface, 
the z = 0 plane. With B * 0, Equation 240 can be simplified, as follows: 

1* . H* ■ 0 


0 » « 0 

A* V 1 

*1 

B* * 0 

O' rot 


(247) 


In this case, only the inversion of 3 x 3 matrices is Involved and this is certainly manageable 
by hand computations. 

A laminated anisotropic composite is governed by 18 independent constants. This number 
can be reduced if symmetry in the method of lamination and symmetry in the constituent 
layers exist. For a homogeneous isotropic plate, the number of independent constants reduces 
to 2. It is important to know that material properties should be referred to the components of 
A, B, and D matrices, or equivalently, the star or prime matrices. They cannot be expressed 
in terms of engineering constants, as homogeneous orthotropic or transversely isotropic 
materials are expressed. 

Because of the complicated coupling effects, due to A^, A 2g , B,, etc., the behavior of a 

laminated composite can best be described using the A, B, and D matrices, or their equivalents, 
without specific references to engineering constants. For example, the effective in-plane 


r 




r 
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shear modulus of a laminated anisotropic material may be Ag g , or 1/Ag f , In general, 


their numerical values are different. Which component is being measured depends on the 
loading condition of the test. For a panel shear test, where only N g jt n, and Nj = = Mj = 0, 



is being measured. In a circular tube under torsion, the loading conditions are as follows: 


Ng t 0, and Nj = N 2 = kj = 0; then A| g is being measured. Component A g6 can be measured 


directly if the loading conditions satis'/ e. / 0, and e. = e„ = k. = 0. These conditions may 
be difficult to achieve. ' 


; ; In the case of uniaxial tension of a laminated composite, say, 0, A^ and a' 12 are the 

components that govern the axial and transverse strains. If a circular tube is loaded along 
' its generator, and A* 2 are related to the axial (longitudinal) and circumferential strains. 

In general, A‘ n t A*^ and a' jo / A* 2< Components A n and A 12 are difficult to measure 

directly because specified strains rather than stresses must be imposed. This is anologous 
to components C,^ being more difficult to measure than S^. 

To avoid confusion, all properties of a laminated anisotropic composite should be reported 
in terms of the component of A, B, and D matrices. The use of engineering constants should 
be avoided. 
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SECTION V 

STRENGTH OF COMPOSITE MATERIALS 


The strength of composite materials, whether unidirectional or laminated, is considerably 
more complicated than the elastic moduli. No unified treatment comparable to that of the 
elastic behavior is available. * ..v 

Three common strength theories can be readily applied to the composite materials. They 
are the maximum stress, maximum strain, and maximum distortional work theories. 

In the unidirectional composite, which is assumed to be orthotropic and quasi-homogeneous, 
the maximum stress theory is expressed by three inequalities: 


V ** 

< 

X 


°y 

< 

Y 

(248) 

°a 

< 

S 



where a , a , and a are the stress components (a state of plane stress 1 b assumed) relative 

x y 8 f. 

to the material symmetry axes; X = axial strength (along the fibers); Y = transverse strength; 
and S = shear strength. Failure of the composite is induced when one or more of the equalities 
in Equation 248 are satisfied. 


The maximum strain theory can also be expressed in terms of three inequalities: 


e y < v. 


(243) 


•• 5 s # 


where e . e , and e_ are the strain components: X - ultimate axial strain; Y « ultimate 
transverse strain; and S ■ ultimate shear strain. Aocordlng to this theory, failure is induced 

when one or more equalities are satisfied. . \\ e 

\■ 

The distortional work theory, in plane stress, can be expressed by 



XY \ v / \ « / 


(250) 


This theory can be represented by a smooth quadratic surface in the stress space. The maxi¬ 
mum stress or maximum strain theory can be represented by a cube in stress or strain space. 
The principal difference between the distortional work and the maximum stress or strain 
theories lies in the existence Interaction among the anisotropic strengths. The distortional 
work theory contains a high level of interaction, whereas the maximum stress or strain theory 
assumes no interaction. Based on available strength data obtained from glass-epoxy com¬ 
posites, the distortional work theory appears to be more accurate than maximum stress 
or strain theory. 
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The strengthofa laminated composite can be predictedif the elastic moduli and the strengths 
of each constituent layer are known, the strain components at each location can be obtained 
from Equation 241, and the stress, from Equation 242. Once these components are known, they 
can be substituted into appropriate strength theories. 

In a laminated composite, not all layers will fail simultaneously. As one or several of the 
layers have failed, the layers which are still intact may be able to sustain the existing load. 
The shifting of the stress distribution within a laminated composite may cause an abrupt change 
in slope in the stress-strain curve of the composite. This is often referred to as the knee. The 
ultimate strength of the composite is reached where the still intact layers cannot carry the 
existing load. This strength analysis of a laminated composite agrees reasonably well with 
available data obtained from glass-epoxy laminated composite behavior. 

A considerable amount of work is still needed for a basic understanding of the strength of 
composite materials. The theories just described are purely phenomenological, and no specific 
reference to the actual mechanisms of deformation and fracture is made. 
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